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Abstract 

c-j . It is shown that the atom-molecule collision problem in the presence of an external electric field can be 

solved using the total angular momentum representation in the body-fixed coordinated frame, leading to 
a computationally efficient method for ab initio modeling of low-temperature scattering phenomena. Our 
calculations demonstrate rapid convergence of the cross sections for vibrational and Stark relaxation in He- 
CaD collisions with the number of total angular momentum states in the basis set, leading to a 5-100 fold 

"53 1 

increase in computational efficiency over the previously used methods based on the fully uncoupled space- 
fixed representation. These results open up the possibility of carrying out numerically converged quantum 
scattering calculations on a wide array of atom-molecule collisions and chemical reactions in the presence 
of electric fields. 
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I. INTRODUCTION 



Recent experimental and theoretical studies have shown that external electromagnetic fields can 
be used as a powerful tool to manipulate molecular collisions and chemical reactivity at low tem- 
peratures [| ll— llVTI . Examples include resonant control of atom-molecule collisions and chemical 
reactions in ultracold molecular gases lll0l- ll2l 1 1411. electric field control of nascent product state 



distributions 115 



16Q, and off -resonant laser field control of motional degrees of freedom [|9|, 



12b. 



These pioneering studies demonstrate that future progress in the field of cold molecules - in par- 
ticular, the ability to create large, dense, and stable ensembles of chemically diverse molecular 
species - will depend to a large extent on our understanding of their collisional properties Jl-7]. 

Theoretical modeling of molecular collision experiments performed at temperatures below 
1 K requires quantum scattering calculations based on multidimensional potential energy surfaces 
(PESs) of unprecedented accuracy, which generally remain beyond the capabilities of modern ab 
initio methods. A way out of this difficulty is to adjust the interaction PESs based on experimen- 
tal measurements of collision observables such as trap loss rates jsi l9l Il8l~ l26|l . The crucial link 
between intermolecular PESs and laboratory observations is provided by quantum scattering cal- 
culations, which yield collisional properties of molecules exactly for a given PES. Because of the 
need to incorporate symmetry breaking effects arising from the presence of external fields [7], such 
calculations are more challenging than their field-free counterparts. In particular, the total angular 
momentum of the collision pair is no longer conserved in the presence of external fields, invalidat- 
ing the standar d appro aches of molecular collision theory based on the total angular momentum 



representation Q27 



28] 



A theoretical formalism for quantum scattering calculations of molecular collisions in external 
fields was developed by Volpi and Bohn and by Krems and Dalgarno [29I 3(3]. The formalism is 
based on the fully uncoupled space-fixed representation, in which the wavefunction of the collision 
complex is expanded in direct products of rotational basis functions and spherical harmonics de- 



scribing the orbital motion of the collision partners in a space-fixed (SF) coordinate frame II29L 
Several groups have used this representation to study the effects of external electric, magnetic, and 



microwave fields on atom-molecule ||llLll2L 



23 



3ll- l34|l and molecule-molecule Q35|— 13VQ collisions 



These studies have shown that the fully uncoupled SF formalism meets with serious difficulties 



when applied to collision problems characterized by strongly anisotropic interactions 115 



35. 



38] 



More specifically, the interaction anisotropy strongly couples different rotational and partial wave 



basis states, leading to very large systems of coupled-channel equations that are beyond the ca- 
pability of present-day computational resources. As most atom-molecule and molecule-molecule 
interactions are strongly anisotropic, this difficulty has precluded converged calculations on many 
interesting collision systems, including Li + HF <H> LiF + H ||15|1 . Rb + ND 3 [26I], Li + NH 113 811 . 
and NH + NH Q37D- 



We have recently developed an alternative approach to atom-molecule and molecule-molecule 
scattering in a magnetic field based on the total angular momentum representation Hi. The total 
angular momentum of the collision complex is approximately conserved even in the presence of 
external fields; thus, using basis functions with well-defined total angular momentum allows for 
a substantial reduction in the number of scattering channels [13911 . This advantage allowed us to 
obtain numerally converged scattering cross sections for strongly anisotropic atom-molecule [40] 
and molecule-molecule 114 ill collisions in the presence of a magnetic field. Magnetic fields interact 
with the electron spin of the molecule, which can be weakly coupled to the intermolecular axis and 
often plays a spectator role during the collision. As a result, while an applied magnetic field shifts 
the energies of the colliding molecules and may lead to the appearance of scattering resonances, it 
hardly affects the mechanism of collision-induced energy transfer. In contrast, electric fields break 
the inversion symmetry of the collision problem and alter the selection rules for parity-changing 
transitions, leading to more dramatic changes in collision mechanisms. Examples include electric 
field-induced molecular states [14211 , dipolar resonances |43|], enhancement and suppression of spin 



relaxation in 2 S and 2 II molecules it 1 1 . 



4411 . and stimulated chemical reactions [8|, 



171 



The purpose of this article is to extend the approach developed in Ref. Ii39ll to describe atom- 
molecule collisions in electric fields. In Sec. II, we formulate the collision problem in the total 
angular momentum representation and outline the procedure of evaluating atom-molecule collision 
cross sections. We then apply our formulation to calculate the cross sections for Stark relaxation 
(Sec. IIIA) and vibrational relaxation (Sec. IIIB) in 3 He-CaD collisions in the presence of an 
electric field. Our results agree well with benchmark calculations based on the fully uncoupled 
SF representation, demonstrating the validity and efficiency of our approach. These findings lead 
us to conclude that numerical algorithms based on the total angular momentum representation are 
a powerful way of carrying out quantum scattering calculations in the presence of electric fields. 
Sec. IV presents a brief summary of main results and outlines future research directions opened 
up by this work. 



II. THEORY 



A non-reactive collision of a diatomic molecule (BC) with a structureless atom (A) in the 



presence of a dc electric field is described by the Hamiltonian (in atomic units) [13911 



H 



2 



(1) 



2/iRdR 2 2fiR 2 

where R is the atom-molecule separation vector, r = rf defines the length and the orientation 
of the internuclear axis (BC) in the SF frame, £ is the orbital angular momentum for the colli- 
sion, V(R, r) is the atom-molecule interaction potential, and \i is the A-BC reduced mass, The 
asymptotic Hamiltonian H mo \ describes the rovibrational structure of the diatomic molecule and 
its interaction with an electric field of strength E oriented along the SF quantization axis Z 

1 d 2 



r _|_ 



+ V(r) - Edcos6 r 



(2) 



2mr dr 2 2mr 2 

where j is the rotational angular momentum, d is the permanent electric dipole moment of the 
molecule with mass m, V(r) is the intramolecular interaction potential 1 4511 . and 6 r is the polar 
angle of the internuclear axis (f) in the SF frame lllll.ll2n. 

The orbital angular momentum £ 2 in Eq. <0Q) can be expressed via the total angular momentum 



of the collision complex J in the body-fixed (BF) coordinate frame as Q28 

I 2 — (J — J? = J 2 + f - J+J- ~ J-h ~ ill 



39] 



(3) 



where J± and j± are the BF raising and lowering operators (note that j± satisfy anomalous com- 
mutation relations Ii46ll ). The BF ^-axis coincides with the vector R and the y-axis is perpendicular 
to the collision plane. 



products of BF basis functions Q28 



As in our previous work [39 1, we expand the wave function of the collision complex in direct 



28,391 



J v, j, k 



(4) 



where k is the BF the projection of J and j, and M is the SF projection of J. In Eq. (J4]), 
\JMk) = a/ (2J + l)/8-ir 2 D(j k (£l E ) are the symmetric top eigenfunctions,D(f2 £ ;) aretheWigner 
D-functions, and are the Euler angles which specify the orientation of BF axes in the SF frame. 
The functions \vjk) = r~ 1 x„j(r)v / 2vrY 7 -fc(#, 0) describe the rovibrational motion of the diatomic 
molecule in the BF frame. The rovibrational functions Xvj ( r ) satisfy the Schrodinger equation 

Xvj(r) = e vj Xvj(r) (5) 



1 d 2 , j(j + 1) , T „ 
+ — — + V(r) 



2m dr 2 



2mr 2 
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where e v j is the rovibrational energy of the molecule in the absence of an electric field |47|]. 
The radial expansion coefficients Ff^ jk (R) satisfy a system of coupled-channel (CC) equations 

d 2 o 
- + 2fiE tot Fj vjk {K) = 2/i 

J',v',j',k' 



dR 2 



Fj M vjk (R) = 2fi (JMk\(vjk\V(R,r,e) + ^(J-j) 



+ H as \J'Mk / }\v'fk'}FP v , fk ,(R), 

(6) 

where E tot is the total energy. The matrix elements of the interaction potential and of £ 2 can 
be evaluated as described in Refs. II28L 13911 . In the absence of an electric field, the asymptotic 



Hamiltonian © has only diagonal matrix elements 

(JMk\{vjk\H as \J'M'k')\v'j'k') = 8jji8 M M<8vv'8jj'hk>£vj (E = 0). 



(7) 



In order to evaluate the matrix elements of the molecule-field interaction in the BF basis, we 



transform the Z-component of vector f to the BF frame [46] 



cos 9 r 



4.7T 



1/2 



101 CY, (p r ) 



(y) V2 E d o;(^)^(^ 



(8) 



The expression on the right-hand contains spherical harmonics of BF angles (9, 0) and Wigner 
D-functions of Euler angles (note that 9 is the Jacobi angle between R and r). Making use of 
standard expressions for angular integrals involving three spherical harmonics 114611 . and neglecting 
the r dependence of d (which is a good approximation for low vibrational states and weak electric 
fields 0480 ) we obtain for the molecule- field interaction matrix element 



(JMk\(vjk\-Edcos9 r \J'M'k')\v'j'k') 

J 1 J' 



x(-) 



M+k-k 



M —M 



-^W'^/[(2J+1)(2J' + 1)(2j + 1)(2j" + 1)] 1 / 2 

j 1 f \( J 1 J' \ ( j 1 f 
/ \ k -q -k' J \ -k q k' 



(9) 



This expression shows that the interaction with electric fields couples basis functions of different 
J. It is because of this coupling that the collision problem can no longer be factorized by symmetry 
into smaller J-subproblems [|28ll . It follows from Eq. © that (i) the external field couplings vanish 
unless J — J' = ±1, and (ii) electric fields couple basis functions of different k, leading to a field- 
induced analog of the Coriolis interaction. Unlike the standard Coriolis interaction, however, the 
interaction with external electric fields couples different fc-states in different J-blocks (assuming 
M = 0). 



The standard asymptotic analysis of the radial solutions to CC equations © at large R gives 
the S'-matrix elements and scattering observables. The analysis proceeds in two steps. First, the 
BF wavefunction is transformed to the SF representation using the eigenvectors of the operator t 2 



51 



52J]. Next, the wavefunction is transformed to the basis in which f/ as is diagonal using the 



eigenvectors of the asymptotic Hamiltonian © in the SF representation. The eigenvalues of H a , 
define the scattering channels and threshold energies e 7 in the presence of an electric field. 
Matching the transformed solutions to the asymptotic form I139ll 



F™{R) -)■ SrfSure 



-i(kyR-£ir/2) 



i \ 1/2 



(10) 



yields the S'-matrix elements describing collision-induced transitions between the channels 7 and 
y with wavevectors k'i = 2fi(E m — e 7 ) = 2fiE c , where E c is the collision energy. The integral 



cross sections can be evaluated from the S'-matrix elements as [3 



3 



39] 



7T 



a 



.2 "}2Y1 l 5 «' 5 77' - S. 
T M 1,1' 



M |2 



(ID 



For the He-CaD interaction, we used a three-dimensional ab initio potential energy surface 
developed by Balakrishnan et al. II45L 15311 . which explicitly includes the r dependence of the 
interaction energy. The rovibrational eigenfunctions Xvj (r) were evaluated by solving the one- 
dimensional Schrodinger equation © using a discrete variable representation (DVR) method [54]. 
The matrix elements of the He-CaD interaction in Eq. © were obtained by expanding the PES 



in Legendre polynomia 



analytically to yield Q28 



s with A, 



3] 



12 and evaluating the integrals over spherical harmonics 



{JMk\{vjk\V(R,r,0)\J'Mk')\v'j'k') = 5j, r 5 kk \{2 3 + l)(2j" + l)] 1 / 2 



^{XvMVxiR^lxv'Ar)} 



A=0 



i A f 

-k k' 




(12) 



The radial coefficients V\(R, r) were evaluated using a 24-point Gauss -Legendre quadrature in 9. 
The r integrals were computed with 30 Gauss-Legendre quadrature points in r e [2.5, 5.6] ao. 

The CC equations © were solved using the log-derivative method [|49ll on a grid of R between 
2 and 100 a with a grid step of 0.1 a . The BF basis set used in Stark relaxation calculations 
(Sec. IIIA) included 10 rotational states (j max = 9); the basis set used in vibrational relaxation 
calculations included 10 rotational states in v = and v — 1 vibrational manifolds of CaD (see 
Sec. IIIB). The cross sections for Stark relaxation were converged to <10%. 



For classification purposes, the eigenvalues of the asymptotic Hamiltonian are assigned physi- 
cal quantum numbers appropriate to a polar diatomic molecule in an electric field: v, j, and m (the 
SF projection of j). In this work, we are interested in low-to-moderate field strengths, where the 
interaction with electric field is small compared to the splitting between the ground and the first 
excited rotational levels. We can therefore keep using j to denote the rotational manifold and m 
to distinguish the Stark states within the manifold, even though j is not a good quantum number 
in an electric field. The assignment procedure works as follows. All eigenvalues of the asymptotic 
Hamiltonian which are close in energy to a particular Stark state \vjm) (that is, |e 7 — e V j m \ < A) 
are assigned the quantum numbers v,j,m. The eigenvalues that do not meet this condition are 
excluded from consideration. In this work, we set A = 0.1 cm^ 1 , however, test calculations show 
that the results are not sensitive to the choice of A as long as E c < A. If this condition is not met, 
problems may arise with distinguishing between elastic and inelastic channels (see Sec. MB). 



HI. RESULTS AND DISCUSSION 

In this section, we first consider the eigenstates of the asymptotic Hamiltonian that define the 
scattering channels in the presence of an electric field (Sec. IIIA). In order to test the performance 
of our approach, we compare the cross sections calculated using the BF total angular momentum 
representation with benchmark calculations based on the fully uncoupled SF representation (Sees. 
IIIB and C). 



A. Asymptotic states 

Figure CD shows the eigenvalues of the asymptotic Hamiltonian © for the ground vibrational 
state of CaD as functions of the applied electric field. The number of total J-states is given by 
Nj = J m ax + 1, where J max is the largest value of J included in the basis set. The eigenvalues 
obtained for J max = 2 and 5 are shown in the upper and lower panels, respectively. The results 
clearly show that H as expressed in the total angular momentum basis has eigenvalues that do not 
correspond to the physical Stark states of the diatomic molecule. This situation is similar to that 



encountered in the case of magnetic fields, and following the terminology introduced in 113911 . we 
will refer to these states as "unphysical". From Fig. Q3 we observe that the number of unphysical 
Stark states increases with the number of J-blocks in the basis set. In addition, the energies of the 

7 



unphysical states become closer to the true Stark energies as J max increases. 

As pointed out before [39j], the origin of the unphysical states shown in Fig. [Dean be attributed 
to the basis set truncation procedure. The total J basis is truncated by restricting the number of 
J-blocks (Nj = J max + 1). However, as follows from Eq. ©, electric fields couple basis states 
in block J to those in block J + 1. When the Hamiltonian matrix is truncated, these couplings 
are left out, resulting in the appearance of unphysical eigenvalues and eigenvectors. In Ref. 0911 
it was shown that the eigenvectors of unphysical Zeeman states are dominated by the largest value 
of J included in the basis set. As a result, the presence of unphysical states has no influence on 
low-temperature collisions in magnetic fields 0911 . 

In order to elucidate the properties of unphysical states, we consider the matrix of the asymp- 
totic Hamiltonian © in the BF basis. In the weak- field limit \Ed\/B e <C 1, we can consider only 
the coupling between the ground and the first excited rotational states in the v = manifold (the 
v index will be omitted for the rest of this section). Arranging the \ JMk)\jk) functions in the 
following sequence: |000)|00), |100)|1-1), |100)|10), |100)|11), |000)|10), 1 100) |00), we obtain 
the matrix of the asymptotic Hamiltonian 



(13) 
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(14) 



(15) 



Diagonalization of Hi yields 



1 



Ai, 2 = B e ± x \Bl + -{Ed)\ A 3 ,4 = 2B e 



(16) 



These energies are the same as those of a polar X S molecule in a dc electric field 113, H. The 
eigenvalues of H 2 

\± = B e ±JBl+ l -(Edy. (17) 
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correspond to unphysical Stark states. The eigenvectors of the unphysical states are given by 



-\Ed B eT jBl + \{Edf 

|A±> = 1000) |00) + V —- 1100)|00) (18) 

L>± L>± 

r I 2 

where D\ = (Ed) 2 + B e - ^ B 2 e =f (Ed) 2 . Eq. (0 illustrates that the field-induced mixing 
between different J-states is proportional to the magnitude of the electric field. Thus, we ex- 
pect that the coupling between the different J-blocks will become stronger with increasing field, 
making it necessary to include more J-blocks in the basis set to obtain converged results even at ul- 
tralow collision energies (see Sec. III). By contrast, the eigenvectors of unphysical Zeeman states 
are, to a first approximation, independent of the field strength 113911 . and so are the convergence 
properties of scattering observables. 

Finally, we note that neglecting the electric-field-induced coupling within the H 2 block leads to 
the disappearance of unphysical Stark shifts ([TBI) . This observation suggests a way to eliminate the 
unphysical states from scattering calculations. Preliminary results obtained with a restricted basis 
set (jmax — 1> ^max = 1) indicate that neglecting the off-diagonal elements of H 2 does provide 
accurate results for both the elastic and inelastic He-CaD scattering. It remains to be seen whether 
or not the procedure can be generalized to larger rotational basis sets. 



B. Stark relaxation in He-CaD (v = 0,j = 1, rrij = 0) collisions 



Figure |2] shows the cross sections for Stark relaxation in 3 He-CaD(w = 0, j = l,rrij = 0) 
collisions calculated using the BF total angular momentum representation. The inelastic cross 
sections are summed over all final Stark states of CaD and displayed as functions of collision 
energy for M = 0. At very low collision energies (in the Wigner s-wave limit) the cross sections 
scale as 1/y/Ec ILlLlZD- At higher collision ener gies , the cross sections display broad oscillations 



due to the presence of scattering resonances 11121 . 



At an electric field of 50 kV/cm, the BF results obtained with J max = 5 are in excellent agree- 
ment with the benchmark calculations over the entire range of collision energies from 10~ 4 cm -1 
to 1 cm -1 . The agreement for J max = 4 is also good at Eq > 0.1 cm -1 . The deviations observed 
above this collision energy occur because the number of total angular momentum states in the 
basis is not sufficient to adequately describe scattering resonances in the entrance and/or exit col- 
lision channels. This is analogous to the lack of convergence at high collision energies observed in 
our previous calculations of atom-molecule collisions in magnetic fields Il39ll . The cross sections 



obtained with J max = 3 are off by ~50 % even in the s-wave regime, which indicates that the 
external field coupling between the J = 3 and J = 4 blocks can no longer be neglected. 

In order to test the performance of our algorithm at higher electric fields, we display in the 
lower panel of Fig. [2]the cross sections calculated for E = 150 kV/cm for different values of J max . 
While J max = 4 cross sections display a similar energy dependence as the benchmark results, 
quantitative agreement requires extension of the basis set to J max = 5. We conclude that it is 
necessary to include more J-states in the basis set to achieve convergence at higher electric fields. 

As shown in the previous section, the properties of unphysical Stark states depend on the mag- 
nitude of the electric field. At higher electric fields the scattering wavefunction contains contri- 
butions from higher J-blocks, making it necessary to increase J max to obtain converged results 
even at ultralow collision energies, as illustrated by the results plotted in Figs. [2]and[3l By con- 
trast, converged results for ultracold atom-molecule collisions in magnetic fields can typically be 
obtained with just two J-blocks 113911 - 



C. Vibrational relaxation: He-CaD(w = l,j = 0,m, = 0) collisions 

In Fig. [31 we compare the cross sections for vibrational relaxation in He-CaD(t; = 1, j = 
0, rrij = 0) collisions calculated using the BF approach with benchmark SF calculations. The 
cross sections are summed over all final rotational states of CaD as plotted as functions of col- 
lision energy for different J max . Balakrishnan et al. considered vibrational relaxation in 3 He- 
CaH(t> = l,j = 0) collisions in the absence of external fields and found it necessary to include 
20 rotational states in the v = and v — 1 vibrational manifolds to achieve numerical conver- 



gence 115 3h . The first excited vibrational state of CaD lies 908.3 cm 1 above the ground state, and 



the rotational constant of CaD is 2.16 cm 1 . In order to properly describe quasiresonant energy 



transfer important at low temperatures 115511 . it would thus be necessary to include at least 20 ro- 
tational states of CaD in each vibrational manifold. A fully uncoupled SF basis with t> max = 1, 
jmax = 20, and Z max = 20 contains 12362 channels. In order to avoid solving large numbers of 
CC equations, we opted to use a restricted SF basis set with v max = 1, j max = 9, and Z max = 9 
to generate benchmark results, which should be adequate for testing purposes provided the same 
convergence parameters v max and j max are used in BF and SF calculations. We emphasize, how- 
ever, that these benchmark cross sections are not physically meaningful (e.g., they may not exhibit 



the quasi-resonance behavior characteristic of vibrational relaxation at low temperatures [53, 
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55|]) 



From Fig. [3] we observe that the BF cross sections obtained for a relatively weak electric 
field (E = 50 kV/cm) are in good agreement with benchmark calculations already at J max = 3. 
Table I demonstrates that a J max = 3 calculation includes only 280 scattering channels, while the 
same calculation performed using the fully uncoupled SF representation requires as many as 1380 
channels. The use of the BF total angular momentum representation thus allows us to reduce the 
number of scattering channels by a factor of 4. The computational cost of solving CC equations 
scales as N 3 with the number of scattering channels Il49ll . so the BF total angular momentum 
representation is more than 100-fold more computationally efficient than the fully uncoupled SF 



representation [|29 



30J]. 



At E = 150 kV/cm, quantitatively accurate results are obtained with J max > 5, while J max = 4 
calculations overestimate the benchmark result by a factor of ~3. Comparison of Figs. |2]and|3] 
suggests that vibrational relaxation cross sections converge more slowly with J max than those for 
Stark relaxation. The gain in computational efficiency (~ 10-fold) is therefore not as dramatic as 
observed for E = 50 kV/cm. Note that the BF inelastic cross sections show an unphysical jump 
at a collision energy of ~0.14 cm -1 . This jump occurs because of the ambiguity of the procedure 
used to assign quantum numbers to unphysical states. As pointed out in Sec. II, the eigenvalues of 
the asymptotic Hamiltonian with energies |e 7 — e V j m | < A are assigned physical quantum numbers 
v, j, and m, where we have chosen A = 0.1 cm -1 . While this procedure works well as long as the 
collision energy is small compared to A, collision-induced transitions between unphysical states 
make it difficult to distinguish between elastic and inelastic channels when this condition is not 
met. This technical difficulty can be eliminated by increasing A or switching to an unphysical 
states-free representation (see Sec. II). 



IV. SUMMARY 

We have presented an efficient theoretical approach to solving the atom-molecule collision 
problem in the presence of an electric field. Unlike previous theoretical work based on the fully 
uncoupled space-fixed representation [3, 3(1, our approach makes explicit use of the total angu- 



lar momentum (J) representation in the body-fixed coordinate frame, in which the atom-molecule 
Hamiltonian has a block-diagonal form in the absence of external fields. The different J blocks 
are coupled only by the molecule-field interaction, making it possible to accelerate convergence of 
scattering observables with respect to the maximum number of rotational states and J-blocks in- 

11 



eluded in the basis set. Our method is thus particularly suitable for quantum scattering calculations 
on atom-molecule (and possibly molecule-molecule) collision systems, where different rotational 
states are strongly coupled by the anisotropy of the interaction pot ential. 

As in the case of molecular collisions in magnetic fields Ii39l1 . truncation of the asymptotic 
Hamiltonian matrix leads to the appearance of unphysical Stark shifts. We have analyzed the 
properties of the unphysical states using a simple 6-state model, which shows that the unphysical 
states arise due to the electric field-induced coupling between different rotational states in adjacent 
J-blocks. The eigenvectors of the unphysical states are linear combinations of different rotational 
and J-states with field-dependent mixing coefficients. Because of the admixture of higher J-states, 
which do not contribute to low-temperature collision observables due to centrifugal barriers, the 
unphysical states are expected to play no role in cold atom-molecule collisions. Furthermore, our 
analytical results suggest that, by neglecting certain coupling matrix elements, it may be possible 
to completely eliminate the unphysical Stark states from scattering calculations. 

In order to test the performance of our method, we applied it to calculate the cross sections 
for vibrational and Stark relaxation in He-CaD collisions in the presence of an electric field. The 
results obtained using the BF approach are in good agreement with benchmark calculations based 
on the fully uncoupled SF representation. Most notably, the number of BF channels required to ob- 
tain converged results is smaller by a factor of 1.5 to 4 (depending on J max ) leading to a 5-100 fold 
gain in computational efficiency (see Table I). These improvements open up the possibility of car- 
rying out highly efficient quantum scattering calculations of strongly anisotropic atom-molecule 
collisions in electric fields, which are of g reat current interest as potential candidate systems for 



sympathetic cooling experiments ||24L 
reactions [8, 



pi g reat 

mm 



26L 40J1 or reactants for electric field-controlled chemical 
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Table I. The number of channels in BF basis sets with different J max for M — 0. All the 
basis sets include 2 vibrational and 10 rotational states of CaD( 2 S + ) and 10 partial waves (for the 
SF basis). The ratio (N SF /N BF ) 3 quantifies the computational efficiency gained by using the BF 
approach. The number of channels in the fully uncoupled SF representation N SF = 1340. 



Jmax 


N BF 


(N SF /N BF ) 


3 


280 


109.6 


4 


420 


32.5 


5 


580 


12.3 


6 


756 


5.6 
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Fig. 1. Stark levels of CaD (bold dashed lines) and eigenvalues of the asymptotic Hamiltonian (full lines) 
as functions of the applied electric field. Upper panel: calculation with J max = 2, lower panel: calculation 
with J max = 5. Unphysical states are shown by red lines. Both calculations are for j max = 9. 



17 




Collision energy (cm ) 



Fig. 2. Cross sections for Stark relaxation in He-CaD(v = 0,j = 1, m = 0) collisions as functions of 
collision energy. The curves are labeled by the maximum value of J included in the basis set (J ma x)> see 
text for details. The electric field is 50 kV/cm (upper panel) and 150 kV/cm (lower panel). The calculations 
are performed for the total angular momentum projection M = 0. Circles - benchmark results obtained 
using the fully uncoupled SF representation. 
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Fig. 3. Cross sections for vibrational relaxation in He-CaD(u = l,j = 0, m = 0) collisions summed over 
all final Stark states of CaD as functions of collision energy. The curves are labeled by the maximum value 
of J included in the basis set ( J max )- The electric field is 50 kV/cm (upper panel) and 150 kV/cm (lower 
panel). The calculations are performed for the total angular momentum projection M = 0. 
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